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In  this  study,  we  obtained  the  approximate  /-states  solutions  of  the 
Schrodinger  equation  interacting  with  a  newly  proposed  potential  known 
as  the  Coulomb-Hulthen-Poschl-Teller  (CHPT)  by  employing  the 
parametric  Nikiforov-Uvarov  method.  The  analytical  expressions  for  the 
energy  eigenvalue  as  well  as  the  vibrational  mean  energy,  vibrational 

specific  heat,  vibrational  free  energy,  and  vibrational  entropy  for 
states  of  carbon  monoxide  (CO)  molecules  have  been  presented  in  detail. 


1.  Introduction 

In  quantum  mechanics,  it  is  possible  to  describe  the  behavior  of  a  physical  system  by  a 
wavefunction  that  can  be  solved  exactly  for  a  given  potential  [1-5],  This  is  because  the 
wavefunction  provides  all  the  necessary  information  about  the  quantum  system.  Knowing  the 
exact  solution  to  the  Schrodinger  equation  plays  a  paramount  role  in  quantum  systems,  but  this 
exact  solution  is  only  possible  for  some  potentials;  thus  approximation  methods  are  usually 
adopted  to  arrive  at  the  solution  [6],  Some  of  the  methods  developed  to  find  the  approximate 
solution  to  the  Schrodinger  equation  include;  the  asymptotic  iteration  method  (AIM)  [7-13],  the 
super  symmetric  shape  invariance  method  [14-19],  the  Nikiforov-Uvarov  (NU)  method  [20-28], 
the  variational  method  [29].  However,  a  good  knowledge  of  these  potentials  helps  in  different 
fields  of  studies  in  physics  such  as  atomic  physics,  condensed  matter  physics,  high  energy  physics 
and  particle  physics  [30-36].  The  solutions  obtained  from  the  combination  of  two  or  more 
potentials  give  significant  results  having  diverse  application  in  quantum  mechanics.  In  one  of  our 
recent  publications  [37],  we  obtained  the  approximate  analytical  solutions  of  the  radial 
Schrodinger  equation  interacting  with  the  Hellmann-generalized  Morse  potential,  which  is  a 
combination  of  the  Hellmann  potential  and  the  generalized  Morse  potential.  Onate  and  Ojonuba 
[38]  while  combining  the  Coulomb  potential,  Yukawa  potential  and  the  inversely  quadratic 
potential  obtained  the  solutions  of  the  radial  Schrodinger  equation.  Their  result  could  be  applied  to 
areas  of  physics  such  as  plasma  physics,  solid  state  physics  and  atomic  physics.  The  analytical 
solution  of  the  Dirac  equation  with  the  Hellmann-Frost-Musulin  potential,  which  is  a  combination 
of  the  Hellmann  potential  and  the  Frost-Musulin  potential,  was  obtained  by  Onate  et  al.  [39].  They 
also  obtained  the  solution  of  the  Klein  Gordon  equation  with  the  combined  potential  (combination 
of  the  general  Manning-Rosen  potential,  the  hyperbolical  potential  and  the  Poschl-Teller  potential) 
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[40].  The  Hellmann-Frost-Musulin  potential  can  be  applied  to  areas  in  condensed  matter  physics, 
atomic  and  molecular  physics,  whereas  the  combined  potential  finds  application  in  areas  like  high 
energy  physics,  nuclear  physics,  atomic  and  molecular  physics.  Ita  and  Ikueba  [41]  solved  the 
Schrodinger  equation  for  the  interaction  of  inversely  quadratic  Hellmann  potential  and  the 
inversely  quadratic  potential  for  any  angular  momentum  number,  /.  They  obtained  the  energy 
eigenvalues  and  their  corresponding  eigenfunctions  in  terms  of  the  Laguerre  polynomials.  Ikhdair 
and  Sever  [42]  obtained  the  approximate  solution  of  the  Dirac  equation  with  a  combination  of  the 
reflectionless-type  potential  and  the  Rosen-Morse  potential  including  the  spin-orbit  centrifugal 
term,  using  the  Nikiforov-Uvarov  method.  In  the  same  vein,  Abdalla  et  al.  [43]  derived  analytical 
solutions  of  the  Schrodinger  wave  equation  for  some  q-deformed  potentials  in  terms  of  Huen 
functions.  While  Falaye  et  al.  [44]  in  their  paper,  solved  the  one-dimensional  Schrodinger 
equation  with  a  combination  of  the  Wood-Saxon  potential,  the  Rosen-Morse  potential  and  the 
symmetrical  double  well  potential,  using  the  Nikiforov-Uvarov  method. 

The  stability  of  a  system  depends  on  the  position  of  particles  within  the  systems.  A  system  with 
vibrating  molecules/particles  is  usually  unstable  because  the  properties  of  the  system  such  as 
entropy,  internal  energy,  free  energy  and  mean  free  energy  change  with  respect  to  the  particles' 
vibration.  Hence,  the  study  of  thermodynamic  properties  becomes  necessary.  Dong  et  al.  [45] 
studied  the  hidden  symmetries  and  thermodynamic  properties  for  harmonic  oscillator  plus  an 
inverse  square  potential.  Also,  Dong  and  Cruz-Irisson  [46]  examined  the  energy  spectrum  and  the 
thermodynamic  properties  for  a  modified  Rosen-Morse  potential.  Baria  and  Jani  [47]  in  2012 
studied  some  thermodynamic  properties  of  liquid  Na,  K,  Rb  and  Cs  at  various  temperature 
measurements.  Oyewunmi  et  al.  [48]  critically  examined  the  thermodynamic  properties  and  the 
approximate  solutions  of  the  Schrodinger  equation  with  Shifted  Deng-Fan  potential  model.  Onate 
and  Ojonubah  [49]  studied  the  thermodynamic  properties  of  a  system  under  the  combination  of  the 
generalized  Poschl-Teller  and  hyperbolical  potentials.  In  2016,  Onate  and  Onyeaju  [50] 
investigated  Dirac  particles  in  the  field  of  the  Frost-Musulin  diatomic  potential  and  the 
thermodynamic  properties. 

In  this  paper,  we  propose  a  new  potential  model  referred  to  as  the  Coulomb-Hulthen-Poschl- 
Teller  (CHPT)  potential  model  which  is  lacking  in  literature.  The  CHPT  potential  model  is  a 
combination  of  the  Coulomb  potential,  the  Hulthen  potential  and  the  Poschl-Teller  potential.  The 
combination  of  these  potentials  produces  a  potential  well,  which  enables  the  molecules  to  be 
trapped  thus  aiding  the  calculation  of  the  vibrational  frequency  of  the  molecules  and  also  the 
thermodynamic  properties  of  the  vibrating  system.  Also,  it  is  known  that  molecule  is 

photo  dissociated  by  absorption  of  vacuum  ultraviolet  photons  on  discrete  lines.  For  this  reason, 
the  thermodynamic  properties  of  states  of  CO  will  be  investigated  using  the  CHPT  potential 
model.  The  CHPT  potential  model  is  given  as: 

X  A  4Ce-"“'^  (1) 

V{r)  = - - 

^  ^  r  l-g-'"' 

where  A,  B,  C  are  potential  strengths,  a  is  the  screening  parameter  and 
t]  =  2a. 

We  employ  the  approximation  scheme  in  Eq.  (2)  suggested  by  Aldrich  [51]  and  Wei  and  Dong 
[52]  in  order  to  obtain  the  I  -  states  solution  of  the  radial  Schrodinger  equation  with  potential  (1) 
and  to  take  care  of  the  centrifugal  term. 
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(2) 


(l-g-'"')' 


2.  The  parametric  Nikiforov-Uvarov  (NU)  method 

The  Nikiforov-Uvarov  (NU)  method  is  based  upon  reducing  the  second-order  linear  differential 
equation  to  a  hyper-geometric  type  equation  [53].  By  introducing  an  appropriate  transformation 
=  ^(x) ,  we  can  then  write  an  equation  of  the  form 


(3) 


■0) 


where  cr(^)  and  cr(s)  are  polynomials  of  degree  two  at  most  and  T  (.s)  is  a  polynomial  of  degree 

one  at  most.  To  use  the  parametric  NU  method,  Tezcan  and  Sever  [54-56]  transformed  Eq.  (3)  into 
the  following  form 


^  \  -A2S  ^ 

!(1-A,s) 


4  2,  3 


Jl-V) 


(4) 


From  the  parametric  NU  method,  the  bound  state  energy  condition  is  obtained  as  [57,  58]  (5) 

-F p  -F -F -F  A^^J A^  ^ -F n{n  —  l) "F  A^  -F  2,^,^  -F  'Z.^jA^Ag  =  0, 

with  the  wave  function  as  (6) 


,  <^3  [  A^q—1, — — /I3Q— 1 1 

'■  4i-2V). 

Where  the  parameters  in  Eq.  (5)  and  Eq.  (6)  are  given  as  follows 


(7) 


2.  =  =  2; +4-„ 

^  =  '^(^  +'^A)  +  ^6’'^o  =  A2  +  2A^-2^,A^^  =  A2-2X^  +  2{^.yjA^ -A^^y 
Al  =  A  ~^/A’A3  =  A  “(VA  “A^/A) 


3.  Radial  solution  to  the  Schrodinger  equation 

The  radial  Schrodinger  equation  is  of  the  form  [38]. 


d^R{r) 


dr^ 


+ 


\E-V{r)- 


Lr 

2// A 


^(r)  =  0, 


(8) 


where  L  =  /(/-f1),  //  is  the  particle  mass,  E  is  the  non-relativistic  energy  and  V(r)  is  the 

interacting  potential.  If  we  substitute  our  proposed  potential  model  in  Eq.  (1)  and  the 
approximation  scheme  in  Eq.  (2)  into  the  Schrodinger  equation  (8),  we  obtain 
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IjuE  Iju 


A  Be 
—  + 
r 


-2ar 


4Ce 


-2ar 


l_e-2- 


(l 


«('■)■ 


(9) 


In  applying  the  NU  method,  a  transformation  variable  of  the  form  s  =  e^°‘’'  is  defined  and 
substituted  into  Eq.  (9)  to  get 


(1-^)  d  1 

■  +  -^ - V  —  +  - 


ds^  5(1-5)  5^(l-5^) 

where 


+  K^S  - 


R[s)  =  Q, 


Kx=- 


2a^n^  ’  " 


=■ 


2a^n^ 


( 2E  +  2  Act  +  5  —  4C ) ,  —  ■ 


2a^h^ 


E  +  2Aa ■ 


V 


2a^n^L^ 

J 


(10) 

(11) 


From  the  comparison  of  Eq.  (4)  and  Eq.  (10),  we  obtain  the  following  parametric  constants 

^2  2^  1, 4.4  0, 2^  — ,  2^  —  +  /q ,  2j  /C2 , 2^  , 

2)  =  —X2  +  /C3  +  /fj  +  —  ,  /l|Q  =  1  +  =  2^1+  +  yj^^. 


(12) 


/l[2  —  .y/^,  /l|3  —  ^  yj^+yl/^. 

Substituting  the  appropriate  parametric  constants  obtained  in  Eq.  (12)  into  Eq.  (5),  we  obtain  the 
energy  equation  of  our  proposed  Coulomb-Hulthen-Poschl-Teller  potential  as: 


E„j  = 


M 


L- 


Afj. 

atd 


2m 


C  + 


a  A  B  a^f2{2l  +  lf  a^t2n{n  +  l) 


2  4 


4m 


2m 


—  n  + 


2jV 


a  ti 


l  +  2n  +  ^j(2l  +1)^  — 

a  H 


(13) 


4.  The  thermodynamic  properties  of  the  Coulomh-Hulthen-Poschl-Teller  potential  (CHPT) 

In  this  section,  we  study  the  thermodynamic  properties  of  CO,  a  diatomic  molecule  with  the 
CHPT  potential  model;  the  energy  of  the  system  is  obtained  as 


E„  = 


2a^h 


2*2 


ju 


^  A,  n  +  T^ 
^  + - 


2(n  +  t) 


J 


(14) 


n  =  0,1,2, . <  [m], 

where 


A.=-^.A,=  ^ 


ah 


\  B  3a^h^  a^h^^ 

Aa - +  6C - 

V  2  SjU  4ju  J 


T  =  ■ 


1  +  jl  “ 


8juC 

24-2 


ah 


(15) 
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We  have  introduced  in  Eq.  (14)  [m]  to  represent  the  largest  inferior  given  by 
m  =  -T  +  ±  for  mathematical  simplicity. 

The  vibrational  partition  function  for  the  CO  molecule  is  obtained  as 


«=0 

where  k  is  the  Boltzmann  constant.  On  substituting  Eq.  (13)  into  Eq.  (16)  yields 

Z(j3,m)  =  ^e  ^ 

n=0 

In  the  classical  limit,  the  summation  in  Eq.  (17)  can  be  replaced  by  an  integral 

=  je  dp,  p={n  +  T). 


Z(P.m) 


4^ 


2u7e ' 


2^A^f3  ^/X^Erfi 


V  ^  7 


4^ 


—  2^1^ 


(16) 

(17) 

(18) 

(19) 


where  (20) 

*  M  A  Al 

Aj  —  2Aj)^,  A4  — 

From  the  vibrational  partition  function  of  Eq.  (19),  the  thermodynamic  properties  of  molecular  CO 
can  now  be  determined  as  follows: 


4.1:  The  vibrational  energy  U: 


U(/^)  =  -—lnZ^M 

op 


r  r 


=  A, 


A,p^ 


44F4{A,) 


-h 


/A3P^(A)A  £> 


-me  ‘ 


where 


K=^^-24^Bvfi 


44^2 


V  ^  7 


-iVX . 


4.2:  Vibrational  specific  heat  capacity  C: 

C(/?)  =  —u  =  -kp""  —U  =  kpA  — ,  i _ _ 

dr  ^  dp  ^  \eP^^p  ^/^(A5) 


1 


{A,A2)  + 


1 


(/!"A,)(A,) 


-2(A,A,) 


M2 


me ' 


where 


(21) 


(22) 


(23) 


(24) 
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A6  = 


V  2  y 


Ay  -  3/  +  ^  o2 


+ 


M2 


Ms 


27Mr 

M3 


A^g ' 


{PKY 


Pm 


M2 


A„  = 


A3p"g^^^^'(A5)A2  g^^^^'cTg 


+ 


MsP 


43/M^ 


2A 


4.3:  Vibrational  mean  free  energy  F: 


F(/?)  =  -/:nnZ„,(/?)  =  -2ln 


P 


2me 


^  2^A^P  ^/^E^fi 


V  ^  y 


•2^/^ 


4.4:  Vibrational  entropy  S: 

S(A)  =  Hnz, „,(/?)  +  yr^lnZ, „,(/?)  = 

oT 


fcln 


2me ' 


27A^^/^Erfi 


V  ^  y 


•\Ia^P  ^Ja^p 


•2^/^ 


UkTA, 


(25) 


(26) 


5.  Some  Theoretic  quantities  and  the  Coulomb-Hulthen-Poschl-Teller  potential  (CHPT) 

In  this  section,  we  calculate  the  Information  energy,  Tsallis  entropy  and  Renyi  entropy  using  the 
probability  density  obtained  from  the  normalized  radial  wave  function.  To  begin,  we  first  calculate 
the  normalized  wave  function.  Given  the  radial  wave  function  as: 

(27) 


[R.AyA=Nly-‘(l-yf[Ff-‘^{l-2y)\  .y  =  e-'-\ 


which  is  equal  to  the  probability  density  p{y),  we  can  easily  calculate  the  theoretic  quantities 


mentioned  above.  Where  a  =  J- ^ ^  +  i{i  +  l)  and  b  =  J(l  +  2£')^  -  .  Normalizing 


2*2 


2a"h 


the  radial  wave  function,  we  have 
\[KPy)^ dy  =  \,  y  =  g'"“. 

Substitute  the  value  of/?„^(y) ,  we  have 


(28) 

(29) 
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a/'2  1  1+v  2 


l-\-  X 


^l-s^ 


where  we  have  used:  2a  =  u,  b  =  v,  x  =  l-2y  and - =  1  -  -  .  Using  integral  of  the  form 


V  ^  y 


n2  ^  _  2'^^  r(t+  n+  l)r(k+  n+ 1) 


(30) 


n  \tY^{t  -\-  k  -\-  fi  -\- \) 


we  have  the  normalization  constant  as: 


N  =2 


1 


oruIuT 

'^2M  +  v  +  3  +  2n^ 

1  2  J 

2U+v+\ 

7  2  r(M  +  n+i)r 

rv  +  3  +  2n^ 

V  2  y 

5.1:  Information  energy 

Information  energy  is  defined  as 

00 

E{p)  =  4;t|  p(r)dr. 

0 

0 

E{p)  =  4;r I p{y)dy.  y  =  exp(-2o;r). 

1 

E{p)  =  4;^  j  p{z)dz.  z  =  1  -  2y. 


-1 


Substitute  the  value  of  the  probability  density,  we  have 

1+v 


E(p)  =  4;rN;,f 


l-z 


-1 ' 


2  y 


V  ^  7 


.n2 


dz. 


Using  integral  of  the  form 

i-ri-yTri+y^'' 


\  ^  j 


Pl"'(y)  'dy  = 


(0_ 

2ric  +  n  +  l)r{d  +  n  +  l) 
n  !r(C"l"  d  2 n+  l)r (c+  d  vi  ~\~V) 


and  the  normalization  constant  calculated  above,  the  information  energy  is  obtain  as 


\(Xi.5AAauY{u  +  n) 

2  2  r(M  +  n  +  l)r 

,  221.68Mri 

hi  iO  1  — 

^2M  +  4n  +  v  +  3^ 

V  2  y 

{u  +  n) 

^\P)  l+v+2u 

2  2  r(M+u+i)r 

^  2u  +  An  +  v  +  3'' 

V  2  y 

(31) 


(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

(38) 


5.2:  Tsallis  entropy 

The  Tsallis  entropy  is  defined  as: 
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T{q)  = 


q-l 


pirYdr 


.  q 


Substituting  for  the  probability  density,  we  have 


T{q) 


1 


q-\ 


I-Ak] 


-1 


N. 


nl 


\  ^  J 


1+v 

^  r 


\  ^  j 


(4T 


dz 


V  J 

Using  integral  of  the  form  in  Eq.  (36),  the  Tsallis  entropy  is  obtain  as 

f 


T{q)  = 


1  12.568 


q-\  q-\ 


AauY{u  +  n) 


\+v+2u 


2  2  r(a  +  n  +  l)r 


2m  +  4m  +  V  +  3 


T{q)  = 


1  12.568 


q-\  q  —  ^ 


J 


8.8192Mr(M  +  n) 


\+v+2u 


2  2  r(M  +  M+i)r 


2m  +  4m  +  V  +  3 


JJ 


5.3:  Renyi  entropy 

The  Renyi  entropy  is  defined  as 

00 

R{q)  = - log4;T[ pir^dr.  q 

1-9  0 

Substituting  for  the  probability  density,  we  have 


R(^)  =  ^— log4;r| 
9-1  _i 


N. 


1+V 

ri-zX-Vl  +  zA-Y 


Ilf. 


\  ^  J 


\  ^  J 


dz. 


V  J 

Substituting  for  the  normalization  constant  with  the  integral  in  Eq.  (36),  the  Renyi  entropy  is 
obtain  as 


R{q)  = 


1.09931 

1-9 


5-2u-v 

2  2  aur{u  +  n) 


r(M  +  M  +  i)r 


2m  +  4m  +  V  +  3 


R(q)  = 


1.0993  17.6384x2 


1+2m+v 

2 


J 


aur(u  +  m) 


1-9 


r(M  +  M  +  i)r 


2m  +  4m  +  V  +  3 


JJ 


(39) 

(40) 

(41) 


(42) 


(43) 

(44) 

(45) 


(46) 


6.  Discussion  of  Results 

Special  cases: 

When  A  =  0,  the  energy  Equation  (13)  becomes 
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^  2a^n^ 

E„,i= - ^ 

M 


L- 


2// 


2t2 


a  h 


B  a^h\2l  +  \f 

a^h^n  +  l)^ 

0  1 
rdv 

2n  ) 

-n2 


(2/  +  1) 


2  8//C 


2t2 


a  h 


2  SjuC 


l  +  2n+  (2/  +  1) 

V  ah 


If  B  =  0,  the  energy  equation  becomes 

2ju 


^  la^h^ 
Kj= - ^ 

M 


An 

ah^ 


a'h\2l+\y  a%'n{n+\) 
2  An  2n 


1 

n  +  - 
V  2y 


(2/  +  ir- 


2  8//C 


-|2 


2t2 


ah 


2  SnC 


\  +  2n+  {21 +  iy- 

V  ah 


If  C  =  0,  the  energy  equation  turns  to 
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The  vibrational  partition  function  for  the  CHPT  molecular  potential  model  has  been  obtained  at 
the  classical  limit.  The  experimental  data  for  molecule  was  taken  from  the 

spectroscopic  experimental  data  reported  [59].  Accordingly,  we  have  adopted  the  following  values 
A  =  1.128A,  B  =  0.149936,  a  =  2.2048A  \  //  =  1. 146x10"^^  g  and  C  =  90531  cm"*  for  the 
potential  parameters  used  in  this  present  study.  In  Figure  1,  the  partition  function  increases 
monotonically  with  an  increase  in  at  a  high  temperature  range.  Also,  the  molecule  starts 
dissociating  at  energy  of  about  5J.  This  can  be  seen  in  Figure  2  with  the  internal  energy  decreasing 
monotonically  with  an  increase  in  p .  The  specific  heat  capacity  of  the  molecule  dissociated  at 
various  values  of  p  except  at  the  value  near  zero  as  shown  in  Figure  3.  This  shows  that  at  low 
temperature  range,  the  spread  of  CO(^A^S^)is  reduced.  The  vibrational  free  energy  of  the  system 

also  shows  some  sensitivity  to  the  range  of  temperature  with  a  negative  value  (Figure  4)  while 
Figure  5  shows  that  the  entropy  of  the  molecule  increases  linearly  with  an  increase  in  p  An  Table 
1,  we  presented  the  numerical  values  for  the  Coulomb-Hulthen  Poschl-Teller  potential.  The 
energy  eigenvalue  decreases  as  the  quantum  number  increases. 


7.  Conclusion 

We  haveinvestigated  the  thermodynamic  properties  of  the  carbon  monoxide  molecule 

using  the  CHPT  potential  model.  The  eigenvalue  obtained  from  the  newly  proposed  potential 
model,  is  then  applied  to  obtain  the  partition  function.  The  explicit  expressions  for  the 
thermodynamic  properties,  such  as:  vibrational  mean  energy  U,  specific  heat  C,  free  energy  F,  and 
entropy  S  of  the  system.  The  effect  of  the  temperature  parameter  p  on  the  various 

thermodynamics  quantities  for  the  A states  of  CO  have  been  investigated  in  detail.  The 
theoretic  quantities  such  as  information  energy,  Renyi  entropy  and  Tsallis  entropy  are  calculated. 
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Figure  1:  Vibrational  partition  function  for  different  values  of  m  , 
as  P  increased  monotonically  with  Z 
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Figure  2:  Variation  of  internal  energy  U  with  p 
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Figure  3:  Variation  of  the  specific  heat  capacity  with  p 
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Figure  4:  Variation  of  the  vibrational  mean  free  energy  with  p 
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Figure  5:  Variation  of  the  vibrational  entropy  with  p 


Table  1:  Energy  spectra  for  Coulomb-Hulthen  Poschl-Teller  potential  for  2p,  3p  and  3d  for  three 
values  of  a. 


State 

a 

fi  =  A  =  -C  =  l 

5  =  2A  =  -2C  =  2 

25  =  A  =  -2C  =  2 

2p 

0.05 

-1.606059113 

-2.270353834 

-1.646249575 

0.10 

-1.598731728 

-2.232544041 

-1.685811919 

0.15 

-1.546183422 

-2.143605317 

-1.688422695 

0.20 

-1.454652035 

-2.011157803 

-1.661349015 

3p 

0.05 

-1.676400130 

-2.329861950 

-1.717237391 

0.10 

-1.747851812 

-2.365182631 

-1.836523320 

0.15 

-1.782073253 

-2.362655259 

-1.925654008 

0.20 

-1.783770598 

-2.328070513 

-1.989020959 

3d 

0.05 

-1.554926971 

-2.211234137 

-1.595883913 

0.10 

-1.404290670 

-2.008428375 

-1.497046307 

0.15 

-1.136625999 

-1.673084350 

-1.296295388 

0.20 

-0.781328125 

-1.240312500 

-1.025078125 
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